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Abstract. The combination of Hipparcos astrometric data with the spectroscopic data of putati ve extrasolar pla nets seems 
to indicate that a significant fraction of these low-mass companions could be brown or M dwarfs (Han et al. 2001). We show 



that this is due to the adopted reduction procedure, and consequently that the Hipparcos data do not reject the planetary mass 
hypothesis in all but one cases. Additional companions, undetected so far, might also explain the large astrometric residuals of 
some of these stars. 
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1. Introduction 

The first attempt to combine the spectroscopic information con- 
cerning extrasolar planetary candidates w ith the Hipparcos data 
has been done by Perryman et al. (1997 ). The goal was to de- 
rive upper mass limits while the radial velocity analysis gives 
the lower limit. 

Han et al. (2001) lately investigated the mass of 30 extra- 
solar companions with period longer than 10 days using the 
Hipparcos (ESA 1997) data. For 27 systems, the inclination 
they derive is smaller than 20° (in 8 cases, it is even less than 
1 °), and they conclude that half of their sample stars might have 
brown or M dwarf secondary. 

In this paper, we first derive similar inclinations for an ex- 
tended set of systems with no limitation based on the period. 
Such a lar ge percentage of small inclinations is very unlikely 
(Sect. 2.1 ). Instead of expl aining it with an observational bias 
of radial velocity surveys (Han et al. 2001), we show that the 
fitting procedure is responsible for the bias. Another model is 
then used to fit the observations. The comparison between the 
two sets of solutions allows us to conclude that the data seldom 
contain enough information to derive a reliable value of the in- 
clination and of the semi-major axis. Such a poor reliability 
will be illustrated in Sect. \. 



Although we show in Sect. |2.2| that the derived inclinations 
are strongly doubtful, we expect a very small number of small 
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inclinations in the putative extrasolar planets sample, and the 
problem is therefore to distinguish between a true astrometric 
motion and the effect of random noise. Statistical tests are then 
proposed in Sect. || in order to detect the real massive secon- 
daries, not due to the effect of random noise or to the perturba- 
tion of other companions. 

2. Statistical properties of the low-mass 
companions sample 

2. 1. The inclinations 

The meth od adopted for fitting the astrometric data (M azeh 
et al. 1999; Halbwachs et al. 2000b| ; Han et al. 2001| ; Zucker 
& Mazeh 2000) generally consists in fixing K\ (reckoned in 
m/s) and the orbital parameters wi, e, P (reckoned in days), 
and T to their spectroscopic values and to impose 



(1) 



The parallax, w, reckoned in mas is adopted from Hipparcos 
and a a is the semi-major axis obtained by astrometry. The 
X 2 thus becomes a 7-parameter expression which is then min- 
imised. From now on, ic will refer to the inclination derived 
using that procedure. 

If the orbital planes of extrasolar companions are randomly 
oriented, Pr(sin i < x) = 1 — V 1 — x 2 . However, when the 
Hipparcos data of 46 systems]] (Table [j]) are fitted using the 

1 Based on the Catalog of extra-solar planets m aintained by 
J. Schneider ( ^ttp://www.obspm.fr/encycl/catalog.htm1 ), content in 



early December 2000 
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Fig. 1. Distribution of sinic (sini^) plotted as open (filled) 
histograms. The dashed and dotted lines represent the ex- 
perimental cumulative frequencies of ic (Campbell) and it 
(Thiele-Innes) respectively. The thick continuous-line gives the 
theoretical cumulative frequency. 




a (mas) sin (i) 

Fig. 2. Inclinations versus a Q sini (Eq. [[]). The triangles and 
squares denote ic and %t respectively. The thick line gives the 
upper bound on ic based on a a smi, assuming a semi-major 
axis of 0.5 mas. The thin line is a least-square fit to the ic 
observed values. 



above algorithm, the distribution of ic peaks close to (Fig. [j], 
open histogram). In order to explain such a distribution, one can 
either argue that (a) the orbital planes of any unbiased sample 
are not randomly oriented; (b) there is a selection effect in the 
sample detected by radial velocity, pushing the inclinations to- 
wards small values, since the binaries have been excluded; (c) 
the inclinations thus derived are plain wrong. We will no longer 
consider explanation (a): even if there could be a preferential 
inclination (e.g in clusters), that it could be (i.e. towards our 
line of sight) would rather be anthropocentric. 

Figure [l] shows that the number of systems with small in- 
clinations (ic) is way too large with respect to the theoretical 
distribution. For instance, there is only a 0.5% probability of 
getting a system with sini < 0.1 whereas 32 systems over 46 
are characterised by sin ic < 0.1. 

Could a selection effect of the radial velocity surveys ex- 
plain this result, as suggested by Han et al. (2001)? To show 
that the answer is in the negative we may take account of all 
companions to solar-type primaries discovered through radial 
velocity studies and compute the probability of getting small 
inclinations. Counting dwarf stars later than F5 from the spec- 
troscopic binaries after Batten et al. (1989] ) and Duquennoy 
& Mayor (1991) and adding the recent low-mass companions 
gives about 450 stars, the exact number does not matter much 
as will be seen below. On the average, we would expect 2 sys- 
tems with sini < 0.1 in our sample, so that the probability to 
get 32 such objects is about 10~ 25 . Even if we allowed for an 
underestimation by a factor 10 (!) of the inclination, to get 16 
systems with sin(10i) < 0.1 would have a 10 -9 probability. 



There may be unpublished spectroscopic companions, so 
let us consider all stars which can be observed. From the Tycho 
Catalogue, there are 72 162 stars with V < 9 and B — V be- 
tween 0.4 and 1.4, corresponding to F5-K8, and only a small 
fraction of this population has been observed indeed. Assuming 
a 11% frequenc y of short period (P < 1000 days) spectro- 
scopic binaries ( jrtalbwachs et al. 2000a ), one can expect on 
the average 0.4 system with sini < 0.01 in this whole popula- 
tion whereas there are 16 syste ms in our sample (resp. 8 in the 
subsample of |Han et al. (200 1| )) i.e. a 10 -6 (resp. 10~ 2 ) prob- 
ability. Facing all these unrealistic probabilities, it is clear that 
the hypothesis (b) above is com pletely ruled out. 

Arenou & Palasi (200& 3 ) and pourbaix (2001 1 ) favoured the 
explanation (c): the inclinations and their formal errors may 
be plain wrong. The right-hand side of Eq. [j] depends on 
the (know n) spectroscopic elements and the parallax only: as 
shown by Halbwachs et al. (2000b), a a is of the order of the 
astrometric precision (about 0.5 mas for the considered stars), 
sinic is then constrained by the spectroscopic a\ sini and the 
parallax: the smaller they are, the smaller the inclination. So, 
the inclinations thus obtained could have no physical meaning 
and could just be artifacts of the fitting procedure. 

2.2. A test of the inclination 

In order to show that these ic are wrong, let us relax some 
constraints on the orbital parameters, namely lo\ and K\. The 
orbital parameters e, P, and T are set to the spectroscopic val- 
ues whereas the others (advantageously replaced by the Thiele- 
Innes constants A, B, F, and G) are derived from the obser- 
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vations. One advantage of this procedure is that, once e, P, 
and T are fixed, the least-square system of equations is linear 
with respect to the parameters and needs no iteration to con- 
verge. Finally the parameters W\, a a , i and are computed 
from the Thiele-Innes constants. We will denote by it the in- 
clinations thus obtained (squares in Fig. |J). If the observations 
do contain the astrometric signature of the planet, the two ap- 
proaches, i.e. Campbell's and Thiele-Innes' ones, both math- 
ematically correct, should yield consistent results. We recall 
here that what we quote 'Campbell's solution' is the determi- 
nation of the Campbell parameters with the supplementary con- 
straint of ai sin i and u)\ to their spectroscopic value. Without 
this constraint the two quoted approaches are identical. 

The uncertainties of %t were derived using Monte-Carlo 
simulation. Instead of comparing each pairs of inclinations, 
we compare the weighted means of the two sets (weight=cr i _2 , 
and using i in the first quadrant). For the Thiele-Innes and 
Campbell sets we obtain respectively it = 83° ±2.7° and 
f c = 0.02° ± 0.004° (6.6 ± 2 for an unweighted mean). The 
two means are clearly discrepant. 

The distribution of It is also given in Fig. [j]. Although it 
looks closer to the theoretical distribution than ic does, one can 
nevertheless notice an excess of edge-on orbits and a deficit of 
intermediate inclinations. The probability of getting no incli- 
nation below 60° in this sample is 10~ 10 . These inclinations 
are as unlikely as those obtained with the Campbell approach. 
Unfortunately, such a bias towards 90° could also be consistent 
with a selection effect due to the spectroscopic investigations. 

Actually, random noise also yields edge-on Thiele-Innes 
orbits. Although, regardless of e, P and T, fitting such a 
null signal leads to Thiele-Innes' elements normally distributed 
around (unbiased), it(A, B, F, G), given by ( fieintz 1978| ) 



2 i T l ( G -A)* + (B + F)* 
2 v (A + G) 2 + (B — Ff ' 



(2) 



follows a distribution centred in it/ 2. When the standard devia- 
tions of the Thiele-Innes elements are all the same, <Ji w 11.7°. 

However, even a likely It would not mean the whole orbital 
solution is sensible. Indeed, even if Campbell's elements are 
replaced with the Thiele-Innes ones, a a still remains closely 
related to the residuals of the coordinates of the star, the reason 
being that the computed semi-major axis roughly follows a law 
proportional to a x(4) law, since a a oc \J A 2 + B 2 + F 2 + G 2 
(from Eq.|). 

If there was some orbital information in the Hipparcos 
data, one would expect a one to one relation between the as- 
trometric Kit (computed using Eq. [I]) and the spectroscopic 
Ki (Fig. ||). However the mean weighted difference between 
the two is 458 ± 62 m/s. What happens here is that, with a a 
about 0.5 mas and an inclination iT close to 90° the astro- 
metric estimate of K\ is almost the inverse of the spectro- 
scopic a\ sin ixzu. Hence, what pushes the inclination towards 
zero with Campbell's approach when a\ sini gets small is also 
responsible for pushing the Thiele-Innes K\ towards a large 
value. 




0.4 0.6 0.8 

a sin i (mas) 

Fig. 3. Ratio of the K\ estimate based on the Thiele-Innes ele- 
ments to the spectroscopic K\ . 



2.3. The law of the semi-major axis 

Since we cannot rely on the inclination, we now turn to the 
semi-major axis. It is indeed their values together with a small 
spectroscopic a\ sin i which pushed the inclinations towards 
small values in the C ampbell approach. 

Han et al. (2001) selected the large secondary masses using 
the ratio ^- so the question is whether large values of this ratio 
are expected or not. We proceed by simulations since one can- 
not easily recover the true ai distribution eith er by assuming 
a Rayleigh-Rice law (|Halbwachs et al. 2000b ) or with an ap- 
proximation formula ( Han et al. 2001 ). The Rayleigh-Rice law 
was shown by simulations to be globally compatible with the 
observed a\ on a sample of stars, but what we need here should 
be done on a star by star basis. 

For each star, we assume the extreme case that all true 
ai are minimum, corresponding to i — 90°, the other orbital 
and astrometric parameters being those obtained in the solution 
with the real data. We then draw at random Gaussian Hipparcos 
abscissae around their expected value and a Campbell solution 
is computed. Using 1000 such simulations we obtain the em- 
pirical probability for each star that the ^ ratio is larger than 
the observed one (column labelled Pr a in Table [j]). 

As can be seen in Fig. |], this probability is not uniform 
(significance level 3%) which is not surprising since the sim- 
ulations assumed no outlier, no other unmodelled astrometric 
perturbation and i = 90°. The main point however is that 
there is no indication of a significant fraction of the real sample 
with a large — ratio, which would correspond to small Pr a in 
the simulations. There are six orbits with a probability smaller 
than 5%, but this is not completely unexpected since this has 
a 3% probability to occur in a 46-sample. Among these stars, 
the result of p CrB and HD 195019 are however less likely 
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Fig. 4. Probability that simulated 24 (assuming i = 90°) are 
larger than the observed one for the 46-sample. 

(prob=0.01%) to occur simultaneously. This could indicate a 
significant orbit or other astrometric perturbations for these 
stars. 

This exercise shows, in some sense, that a whole sam- 
ple made of edge-on orbits would be more probable than the 
large fraction of face-on orbits which are obtained by the direct 
Campbell approach. Edge-on orbits for all stars are of course 
unrealistic, and this shows merely that there is no orbital infor- 
mation in the Hipparcos data for the sample as a whole. 

3. Screening the astrometric solutions 

At this point, it is clear that there is no indication in the plan- 
etary sample as a whole of a significant fraction of low incli- 
nations/large secondary mass. Due to the increasing number of 
planetary candidates, it is however expected on a statistical ba- 
sis that a bunch of them may have a mass much larger than the 
minimum mass. 

Since there is evidence that a direct use of the orbital pa- 
rameters (without taking their statistical distribution into ac- 
count) may lead to incorrect conclusions, we need a set of tests 
which can tell us whether the Campbell results are correct or 
not. Since by construction the Campbell analysis is consistent 
with the spectroscopic orbit, this is where the Thiele-Innes ap- 
proach reveals useful. 

The tests described below intend to show whether an orbital 
model improves the fit of the Hipparcos data or not, whether 
the Thiele-Innes elements are consistent with Campbell's ones 
or not and, finally, whether the astrometric and spectroscopic 
orbits agree with each other or not. This is an essential point: 
before combining the spectroscopic and astrometric data, as is 
done in the Campbell approach, one has to be sure that both 
solutions are compatible. If all these tests give a positive result, 



then the orbital elements obtained with the Campbell approach 
may be trusted. 

In the other case, one may either conclude that the or- 
bital elements were obtained by chance only, or that there is 
a perturbation in the astrometric solution due e.g. to another 
companion. The presence of other companions around several 
stars hosting planetary candidates can be detected from a ra- 
dial velocity residual trend (Fischer et al. 2001) and some of 
these companions are being discovered by adaptive optics: HD 
1 14762, t Boo (|LLoyd et al. 2000| ) and Gl 86 ( fisetal. 200)\ ). 
Arenou & Palasi (2000) combined the spectroscopic data with 
the Hipparcos and Tycho2 data (H0g et al. 2000) for the plan- 



etary candidates and the use of long term proper motion gave 
a hint of possible other companions (period larger than several 
decades). 

The spectroscopic orbital elements are generally precise. 
So, from now on, we treat them as constants. As we have seen, 
the statistical distributions of the semi-major axis and of the 
inclination are not trivial, so our tests must be based on better 
known parameters. This is the case for the A, B, F, G Thiele- 
Innes parameters: obtained by a linear least-squares solution, 
they may be considered Gaussian around 0, as has been con- 
firmed by direct simulations. Concerning real Hipparcos data, 
it has been shown that the errors on the astrometric p arame- 
ters could be considered Gaussian ( Arenou et al. 1995 ). So the 
normality hypothesis for the Thiele-Innes elements is not just 
a convenient hypothesis but is plainly justified. Of course, a 
departure from normality may be due to outliers. Apart from 
possible instrumental errors, however, outliers may be normal 
observations when a good orbital model is used in place of the 
standard "single star" model adopted for the major part of the 
Hipparcos stars. 

3. 1. Need for an orbital motion 



As pointed out by Pourbaix (2001 ), none of the Campbell-like 
orbital solutions but four does improve the fit of the Hipparcos 
intermediate astrometric data (IAD). Could it be different with 
a Thiele-Innes solution? Do the two additional degrees of free- 
dom significantly reduce the x 2 ? The quantity 



F = 



N-9x% 



Xt 



(3) 



follows a F-distribution with (4, N — 9) degrees of freedom 
( Bevington & Robinson 1992 ). N denotes the number of data 



points and \ s and Xt ^ me va l ue of the \ w i m the 5- 
parameter (single star) model and Thiele-Innes' orbital model 
respectively. We reject the null hypothesis 'no orbital wobble 
present in the IAD' if the probability of getting F larger than 
F(4, N — 9) is lower than 5%. This P-value is given as Pi - 2 in 
Table [[] wh ere we also recal l as Pri the analogous probability 
obtained by Pourbaix (2001 ) for the Campbell approach. 
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Table 1. List of systems ( pourbaix 2001 ) and the statistical tests. Pr Q =Probability to get a a a /a a larger than the observed 
value assuming i=90°; Pri=Pr(F c > F(2,N - 7)) (|Pourbaix 200l|); Pr 2 =Pr(F > F(4,N - 9)); Pr 3 =Pr(xl s > X 2 (4)); 



Pr 4 

\D\ < 



Pr(x!c 



> x 2 (4)); The Campbell approach may be accepted when the conditions Pri < 5%, Pr 2 < 
Pr 4 > 5%, Pr 5 > 5% are fulfilled. 



5%, 



> 
Pr 3 



< 5%, 



HIP 



HD/name Ref 



Pr a 


Pri 


Pr 2 


Pr 3 


D 


Pr 4 


Pr 5 


(%) 


(%) 


(%) 


(%) 




(%) 


(%) 


67 


66 


89 


89 


0.24 


99 


87 


14 


51 


12 


33 


0.77 


44 


15 


55 


62 


52 


50 


-0.38 


66 


30 


16 


23 


14 


14 


0.93 


41 


14 


37 


20 


25 


25 


-0.25 


71 


34 


2 


6 


19 


5 


0.53 


89 


61 


79 


36 


65 


51 


0.24 


96 


73 


39 


75 


64 


64 


0.44 


72 


36 


50 


46 


59 


59 


0.59 


86 


52 


28 


20 


28 


23 


0.34 


95 


73 


42 


40 


60 


61 


0.21 


92 


63 


36 


44 


71 


59 


-0.24 


95 


71 


86 


55 


58 


23 


0.03 


41 


13 


3 


8 


8 


1 


-2.07 


48 


26 


46 


46 


93 





0.38 








54 


59 


77 


74 


0.11 


93 


66 


30 


30 


13 


12 


-0.65 


32 


9 


25 


61 


68 


61 


0.04 


81 


46 


5 


13 


15 


9 


-0.07 


52 


20 


2 


3 


8 


4 


-1.04 


77 


40 


86 


35 


62 


42 


0.60 


93 


69 


26 


22 


31 


27 


-0.17 


74 


38 


97 


86 


98 


96 


-0.12 


99 


90 


56 


41 


52 


57 


0.21 


86 


52 


38 


42 


69 


59 


-0.58 


95 


71 


60 


55 


74 


34 


-0.72 


25 


14 


21 


27 


47 


35 


-0.51 


87 


54 


46 


61 


40 


29 


0.44 


44 


15 


59 


55 


84 


90 


0.26 


99 


90 


.2 


.2 


1 


.1 


-0.34 


99 


87 


36 


25 


37 


44 


0.19 


86 


45 


44 


50 


96 


77 


0.22 


95 


75 


88 


60 


92 


63 


-0.17 


98 


91 


13 


10 


30 


37 


-0.02 


97 


79 


15 


19 


28 


33 


0.69 


81 


46 


44 


41 


65 


67 


-0.52 


95 


70 


34 


17 


46 


62 


-0.22 


99 


95 


12 


9 


24 


22 


-0.27 


94 


67 


30 


73 


48 


20 


-2.48 


30 


8 


.5 











-4.68 


6 


1 


48 


19 


15 


8 


-2.08 


40 


13 


53 


41 


72 


66 


0.00 


99 


97 


49 


49 


43 


53 


0.32 


74 


37 


46 


65 


86 


79 


-0.13 


99 


95 


61 


54 


60 


68 


0.62 


86 


53 


52 


10 


17 


28 


0.32 


82 


47 



1292 

5054 

7513b 

7513c 

7513d 

8159 

9683 

10138 

12048 

12653 

14954 

16537 

26381 

27253 

31246 

33719 

43587 

47007 

47202b 

47202c 

50786 

52409 

53721 

64426 

65721 

67275 

68162 

72339 

74500 

78459 

79248 

89844b 

89844c 

90485 

93746 

96901 

97336 

98714 

99711 

100970 

108859 

109378 

113020 

113357 

113421 

116906 



GJ3021 

6434 

v And b 

v And c 

v And d 

10697 

12661 

G186 

16141 

HR810 

19994 

e Eri 

37124 

38529 

46375 

52265 

55 Cnc 

82943 

83443b 

83443c 

89744 

92788 

47 UMa 

114762 

70 Vir 

t Boo 

121504 

130322 

134987 

pCrB 

14 Her 

168443b 

168443c 

169830 

177830 

16 Cyg B 

187123 

190228 

192263 

195019 

209458 

210277 

G1876 

51 Peg 

217107 

222582 



Naef etal. (2000) 



Queloz et al. (2000a) 



Butler etal. (1999 ) 



Butler etal. (1999 ) 



Butler etal. (19991) 



Vogt et al. (2000) 



Fischer et al. (20011) 
Queloz et al. (2000b|) 



Marcy et al. (2000|) 



Kiirstcr et al. (200CD 



Queloz et al. (2000a) ) 



Hatzes et al. (2000| ) 
Vogt et al. (20001) 



Fischer et al. (2001h 



Marcy et al. (200C ) 



Butler et al. (200C ) 



Butler etal. (1997) 



Naef etal. (2001) 



Mayor et al. (200C ) 



Mayor et al. (2000,) 



Korzennik et al. (2000) 



Queloz etal. (2000a) 



Butler & Marcy (1996| ) 



Vlarcy etal. (1999) 



Marcy & Butler (1996| ) 



Butler etal. (1997) 



Queloz et al. (2000a) 



LJdry etal. (2000a) 



Vogt et al. (200C \ 



Noycsct al. (199S ) 



LJdry ct al. (20001: ) 



Udry et al. (20001: ) 



LJdry et al. (2000b; ) 



Naef etal. (2001 



Vogt et al. (200' 



Cochran et al. (19971 ) 



Vogt et al. (200CD 



Sivan et al. (200C ) 



Vogt et al. (200C ) 



Vogt et al. (200CD 
Vlazeh et al. (2000) ) 



Vogt et al. (200(1" 



Dclfosse etal. (1998) 



Mayor & Queloz (1995p 



Vogt et al. (200C ) 



Vogt et al. (200C ) 



3.2. Significance of A, B, F, andG 

Getting a substantial reduction of the x 2 with the Thiele-Innes 
model does not necessary mean the four constants are signifi- 
cantly different from 0. In order to assess that, we compute 



where X = (A, B, F, G) and V is the covariance matrix of A, 
B, F, and G.'t' denotes the transposition. xls follows a x 2 - 
distribution with 4 degrees of freedom when there is no orbital 
information. When all formal errors on Thiele-Innes parame- 
ters are identical, this xls i s proportional to a a . We then reject 
the null hypothesis 'a a (mas) is significantly different from 0' 



Xls = X'V^X 



(4) 
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if the probability that x\s excee ds X 2 (4) is less than 5%. This 
P-value Pr 3 is given as in column 7 of Table |5J 

3.3. Consistency between spectroscopy and 
astrometry 

Testing such a consistency is a bit tricky. Indeed, the nature of 
the information supplied with by both sides is essentially dif- 
ferent. The astrometry should however be able to recover a sin i 
and oji, both given by spectroscopy. There are a few useful re- 
lations based on the Thiele-Innes elements such as 



L = A 2 + B 2 — (F 2 + G 2 ) = (a a sini) 2 cos(2wi) (5) 



and M = A 



B 



F 2 + G 2 = 2a 2 a - (a a sini) 5 



that can be used to build up a test. For instance, 



D 



Li — L r 



2a\/M r + 2ct 2 



(6) 



(7) 



follows a distribution of null expectation and unit variance. Li 
denotes the left-hand side of the definition of L, i.e. based on 
the Thiele-Innes constants only while L r and M r are computed 
from the spectroscopic elements and the parallax. In the latter, 
the true a a is unfortunately unknown, so it is derived from the 
Campbell-like solution, a stands for the mean standard errors 
of A, B, F, and G. Because we neglect the correlation between 
the Thiele-Innes constants, and with the a a deduced from the 
Campbell analysis, D is not Gaussian, and \D\ < 2 does not 
strictly correspond to a 95% confidence level. Indeed, as can be 
seen in Table |], the distribution of D is leptokurtic: a a being 
overestimated, \D\ is generally underestimated. A large value 
of \D\ is then a good indication that the astro me trie solution is 
not consistent with the spectroscopic orbit. 

3.4. Consistency between Thiele-Innes' and 
Campbell's solutions 

Since both methods are applied to the same data set, they 
should also give similar results. We have already shown that 
edge-on orbits are consistent with pure noise signal when fit- 
ted with the Thiele-Innes parameters while Pourbaix (2001 ) got 
similar results for face-on orbits with Campbell's elements. 

Here also, designing a test may be difficult, especially be- 
cause the distribution of £1 is anything but symmetrical. So, 
instead of comparing the two sets of Campbell elements, we 
compute their Thiele-Innes expressions. We then build 



Xlc = Y'V^Y 



(8) 



where Y = (A - A c , B - B c , F - F c , G - G c ), V has the 
same meaning as in xls (Eq. |j) and (Ac, Be, Fc,Gc) are 
computed from the Campbell elements. If we neglect the cor- 
relations between Campbell and Thiele-Innes, the quantity xlc 
follows a x 2 -distribution with 4 degrees of freedom. We reject 
the null hypothesis 'Thiele-Innes' and Campbell's approaches 
yield consistent solutions' if the probability that xlc excee ds 
X 2 (4) is less than 5%. That probability Pr4 is given in Table 



[I]. Because of the correlations between the A, B, F, G obtained 
by both methods, a rejection at a 5% level will occur in much 
less than 5% of the cases. In fact, xlc i s probably closer to a 
X 2 -distribution with 2 degrees of freedom, so this test is less 
discriminant than the next one. 

3.5. A small inclination 

Finally, one may answer both whether a small inclination is 
present or whether the astrometric solution is compatible with 
the spectroscopic one, still using the Thiele-Innes elements as 
an external control. Specifically, 



A — G 
B + F 
and 



wa\ sini cos(f2 — u>i) tan( 
wai sini sin(f2 — u>x) tan( 



(9) 



A + G = max sini cos(ri + wi) cot(-) 



B-F = 



waisini sin(fi + Wi) cot( — ) 



(10) 



In the case of the planetary candidates, a sini is small, 
and, if the inclination is really small, then the right-hand sides 
of Eqs. | are null; when the inclination is found to be near 
180°, the Eqs. |k] should be « 0. We test whether Eqs. g (if 
cosic > 0) or Eqs. 10 (if cosip < 0) are verified. If Z de- 



notes this two-dimensional vector, and W its covariance matrix 
computed using the covariance matrix of (A, B, F, G), then 



x\ = z'w-'z 



(11) 



should follow a x 2 with two degrees of freedom under the 
null hypothesis of a Thiele-Innes solution compatible with the 
Campbell/spectroscopic solution. The P-value Pr 5 is given in 
Table |l|. We will reject the compatibility between Thiele-Innes 
and Campbell if Prs is smaller than 5%. 

3.6. Discussion 

Since the proposed tests are not independent, we cannot com- 
pute an overall probability. While the first three tests check 
whether the astrometric solution is firmly established, the last 
three ones (useless in the case of a 'noise only' orbit) assess the 
agreement between the astrometric and spectroscopic orbits. In 
fact Pr! and Pr 5 would probably be enough for investigations 
of future planetary candidates, provided that a signal is indeed 
present at the spectroscopic period. It appears that, with the ex- 
ception of p CrB, almost all stars in the present sample fail to 
at least two of the proposed tests. 

Under the assumption of a 'noise only' orbit, one expects 
that the computed probabilities are uniform on this sample. 
Once again, a departure of the Hipparcos residuals from strict 
normality (e.g. due to outliers or long period companions) 
could produce spurious solutions which in turn would contam- 
inate the probabilities. Nevertheless, using a Kolmogorov test, 
the uniform null hypothesis is accepted for Pr 2 (P-value=88%), 
Pr 3 (7%), Pr 5 (98%). 
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If this was necessary, this terminates the demonstration that 
the Hipparcos data mainly contains noise. This does not mean 
that there is no brown or red stars unfortunately interpreted as 
planet (apart p CrB), but rather that Hipparcos is of little help to 
demonstrate it, to the contrary of the Han et al. (2001 ) analysis. 



4. Individual stars 

4. 1. Validity ofi c with HD 209458 

HD 209458 (HIP 108859) is, so far, the only case where the 
inclination is accurately known. Indeed, transits of the planet 
have been monitored since 1999 (Charbonneau et al. 2000), 



thus setting a strong constraint on the inclination: i = 86.1 ± 
1.6° dMazeh et al. 20001). 



Hipparcos observed that system and it has been noticed 
later on that it did detect the planet . . . thanks to its photo- 



metric signature (Soderhjelm 1999; Robichon & Arenou 2000; 
Castellano et al. 2000| ). However, that result does not mean the 
astrometric signature is also in the IAD. Thus, when ic and 
£1 are fitted together with the five astrometric parameters, the 
inclination is i c = 0.019 ± 0.0097°. Fitting the Thiele-Innes 
elements yields an inclination %t = 76 ± 18°. 

It appears that the 'noise only' null hypothesis can't be re- 
jected as indicated by three tests (Pri, Pr2 and P1-3). Moreover 
the hypothesis of a small inclination has a small P-value 
(Pi - 5=13%). On this real case where the true orbit is known, 
the proposed tests thus seems discriminating enough to detect 
wrong Campbell solutions. 

42. p CrB 

The suggesti on that the secondary c ould in fact be an M dwarf 
was m ade by Gatewood et al. (2001 ) using both Hipparcos and 
MAP ( Gatewood 1987 ) data, and the solution in the Hipparcos 
Catalogue was already detected as orbital with a period twice 
the spectroscopic period. As shown by Gatewood et al. (2001) 
a periodogram exhibits a strong signal at this period. 

This is the only star in the 'planetary' sample which suc- 
ceeds to all tests: while the hypothesis of a 'noise only' orbit 
is significantly rejected by the three first tests, there is no indi- 
cation that the spectroscopic, Campbell and Thiele-Innes orbits 
could be inconsistent, using the other tests. Consequently, the 
main result of Gatewood et al. (2001) may not be rejected. 




Period (d) 



Fig. 5. Periodogram of HD 10697 assuming e and T from the 
known companion. 



be more likely, but would have been seen in the spectroscopic 
data. 

Repeating this periodogram with the Campbell analysis 
is shown in Table [| While the obtained semi-major axis in- 
creases with trial period, the inclination decreases accordingly, 
the same large signal to noise w 3 is obtained, and the goodness 
of fit (GOF) is not improved at the spectroscopic period. The 
Campbell analysis simply adjusts itself to the constrained spec- 
troscopic elements, without showing any sensible information 
in the astrometric data. 

This explains the negative results of our tests in Table |l[ 
Due to the size of the astrometric residuals (partly independent 
from the orbital perturbation), the obtained semi-major axis is 
probably overestimated, but it is difficult to know to which ex- 
tent. Consequently the exact mass of the secondary cannot be 
obtained at face value using a combination of the Hipparcos 
data with the spectroscopic orbital elements. 



4.3. HD 10697 

Since the minimum reflex semi-major axis is about 0.4 mas, 
and thus suggest s that Hipparcos could g ive a hint of the true 
semi-major axis, Zucker & Mazeh (200C ) combined the astro- 
metric and spectroscopic data and obtained a a = 2.1 ±0.7 mas, 
corresponding to a 38 ± 13 Mj mass for the secondary. This 
result would indicate that the secondary is in the brown dwarf 
domain. 

Although the signal should be present to some extent in the 
Hipparcos data, a periodogram based on the unit-weight error 
of the Thiele-Innes solution shows no signal at the 1072 days 
spectroscopic period (Fig. |]). Instead a 345 days period would 



Table 2. Campbell analysis with different trial periods on 
HD 10697 



Period 


a a 


GOF 


(days) 


(mas) 




700 


1.9 ±0.7 


1.31 


800 


1.9 ±0.7 


1.28 


900 


1.9 ±0.7 


1.27 


1000 


2. ±0.7 


1.32 


1072 


2.1 ±0.7 


1.28 


1100 


2.2 ±0.8 


1.25 


1200 


2.7 ±0.9 


1.24 


1300 


3± 1 


1.29 
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Fig. 6. Periodogram of HD 195019 assuming e and T from the 
known companion. 



satellite scan orientation) which could be wrongly attributed to 
the reflex motion due to the low-mass companion, especially if 
the companion had a significant motion during the Hipparcos 
observation period. 

4.5. HD 38529 

A periodogram reveals nothing special at the 14-day spectro- 
scopic period. However, this is another star for which Han et al. 
(2001) suggest that the secondary could be an M dwarf. They 
also quote that this star received an acceleration solution in the 
Hipparcos Catalogue. 

Before looking to the small signal due to the spectroscopic 
companion, one has to remove the effect due to the long pe- 
riod companion. When the acceleration terms in time are thus 
subtracted, the semi-major axis of the Campbell solution be- 
comes 1.3 ± .7 mas, i.e. not significant, whereas it was 2.1 ± .6 
mas before; when a cubic instead of quadratic term in time is 
accounted for, a a — .9 ± .7 mas. 

The influence of a long-period companion, also detected 
in the spectroscopic data (Fischer et al. 2001 ), is thus a much 
more natural explanation to the size of the Hipparcos residuals 
than a tiny inclination for the short period companion. 



4.4. HD 195019 

This is probably a typical case of how one can be mislead by 
a direct Campbell analysis: the fit of the IAD is significantly 
improved by the orbital model but fails the last two tests, i.e. 
the astrometric orbit is clearly not consistent with the spectro- 
scopic one. The periodogram based on the orbital parameters 
of the known companion (Fig. ||) reveals a peak at 68 days 
and nothing around 18 days (the orbital period of the planetary 
companion). 

We cannot find any satisfactory solution by fitting the 12 
parameters with a global optimisation technique. We end up 
with a minimum \ 2 which corresponds to e « 1 and a total 
correlation between F, G, and e. Even if that solution does 
minimise x 2 (12), it is physically very unlikely. 

Instead of fitting the 12 parameters in once, one can use the 
periodogram and build up the solution by adding the parame- 
ters one by one. Thus, by first fixing P to 67.935d, one derives 
T = JD 2, 448, 296.8 which yields the lowest value of the 
unit-weight error. The resulting orbit would have an inclination 
of 65± 18° and K i = 24±8km/s. Unlike the orbit of the known 
companion, this solution successfully passes the last two tests. 
However such a solution would have been present in the radial 
velocity residuals. 

Another explanation could originate in the fact that HD 
195019 is known since 1881 as a binary star, with a B com- 
ponent 4.5" apart (3.5" in 1988) with a position angle of 
323°(resp. 330°), from the WDS Catalog. Despite this fact, 
it has been dealt as a single star in the Hipparcos Catalogue 
(although a 'suspected non-single star' flag is set), because no 
convincing double star solution could be used. Although w 3 
mag fainter, the influence of the secondary may have produced 
astrometric residuals (the projected separation varying with the 



4.6. HD 83443 

There is also no signal around the periods of the two planetary 
candidates (w 3 and 30 days) in the Hipparcos data. Instead 
a period of 6, 79 or 462 days (Fig. [^) would be much more 
likely, but it would clearly have been seen in the spectroscopic 
data. An unmodelled astrometric perturbation due to a longer 
period companion would provide a more logical explanation, 
since there is a possible trend in the radial velocity residuals 
(S. Udry, private communication). 

Still, based on the tests we presented, the astrometric solu- 
tion for HD 83443c was about to be accepted. The reason for 
this is the large a a = 2 ± 0.7 mas. 

One may ask whether any other period could have given 
the same signal as the one obtained with a 30 days period. All 
above mentioned tests were based on normality assumptions 
about the astrometric residuals but another, non-parametrical 
test may also be used. A periodogram, based on the amplitude 
of the signal (root-sum-square of the 4 Thiele-Innes parame- 
ters) has been done: using 5000 trial periods between 2 and 
1500 days, equally spaced in frequency, a Thiele-Innes solu- 
tion is performed. It happens that 30% of the trial periods give 
an amplitude larger than the one obtained at the period of the 
spectroscopic companion. In other words, the astrometric result 
has no special meaning. 

5. Conclusion 

In summary, the available Hipparcos data is not precise enough 
to show that a significant part of the planetary candidates could 
be of stellar nature. 

Although we expect a very small number of small incli- 
nations in the putative extrasolar planets sample, the inclina- 
tions derived using the Campbell elements together with the 
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Fig. 7. Periodogram of HD 83443 assuming e and T from the 
second companion. 



Hipparcos data are pushed towards unreasonably small values. 
As we have shown, this is not due to an observational selection 
but rather to the employed reduction method and probably also 
due to the presence of longer period, unresolved companions. 

Though the Thiele-Innes alternative method gives inclina- 
tions more in line with what is generally expected from the 
spectroscopic detection, the estimated astrometric semi-major 
axis still does not provide a statistically unbiased estimate of 
the true reflex motion. Indeed, in the case of pure noise in the 
astrometric data, the Campbell approach would give inclina- 
tions biased towards 0° while Thiele-Innes recovers inclina- 
tions biased towards 90°. So, none of the two methods can be 
trusted, except when it can be shown that they agree together, 
which is the basis of one of the tests we propose. The other tests 
determine whether the astrometric orbit is significant or might 
result from noise only. 

It must be pointed out that the bias on the semi-major 
axis and inclination only occurs for negligible orbital wobbles. 
When a true orbital information is present in the Hipparcos as- 
trometric data, then the orbital parameters are much better be- 
haved. For instance, we may question whether the rejection of 
putative brown dwarfs into the stellar domain by the Halbwachs 
etal. (2000b) study was justified. The stars HIP 13769, 19832, 
62145, 63366 and 1 13718 which were found with a mass above 
the H-burning limit with a 2a significance and HIP 77152, one 
a above, have been tested using the tests described above. All 
these stars succeed to all tests, apart from HIP 19832 which 
fails to the last test (significance=3%, this may be by chance 
only since the reality of the orbit is clear, as seen e.g. by pe - 
riodogram). The main conclusion of [Halbwachs et al. (2000b ), 
the deficit of short period brown dwarf secondaries, is thus con- 
firmed by this analysis. 



In contrast, in the sample of putative planets, there is no 
indication that a significant fraction of secondaries could be 
brown or red dwarfs, with the possible exception of p CrB. Is 
that to say that the Hipparcos astrometry is useless for the ex- 
trasolar planet study? Not completely, since at least an upper 
limit on the secondary mass may be obtained, and Hipparcos 
may also give a hint of other, longer period companions. For a 
satellite which was originally planned to get positions, proper 
motions and parallaxes, this is not negligible. 
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